Simulating field-aligned diffusion of a cosmic ray gas

The macroscopic behaviour of cosmic rays in turbulent magnetic fields is discussed. An implementation of anisotropic diffusion of cosmic rays with respect to the magnetic field in a non-conservative, high-order, finite-difference magnetohydrodynamic code is discussed. It is shown that the standard implementation fails near singular X-points of the magnetic field, which are common if the field is random. A modification to the diffusion model for cosmic rays is described and the resulting telegraph equation (implemented by solving a dynamic equation for the diffusive flux of cosmic rays) is used; it is argued that this modification may better describe the physics of cosmic ray diffusion. The present model reproduces several processes important for the propagation and local confinement of cosmic rays, including spreading perpendicular to the local large-scale magnetic field, controlled by the random-to-total magnetic field ratio, and the balance between cosmic ray pressure and magnetic tension. Cosmic ray diffusion is discussed in the context of a random magnetic field produced by turbulent dynamo action. It is argued that energy equipartition between cosmic rays and other constituents of the interstellar medium does not necessarily imply that cosmic rays play a significant role in the